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Abstract 

We present a perturbation theory to find the response of an anisotropic 
DNA to the external tension. It is shown that the anisotropy has a 
nonzero but small contribution to the force-extension curve of the 
DNA. Thus an anisotropic DNA behaves like an isotropic one with an 
effective bending constant equal to the harmonic average of its soft 
and hard bending constants. 

1 Introduction 

One of the most successful theories to describe the physical behavior of a 
long DNA molecule is the elastic rod model [1]. In this theory, the DNA is 
modeled as a continuous rod with intrinsic twist (to account for the helical 
structure of DNA) which changes its conformation in response to external 
forces or torques. The response of the DNA to an external stress is then 
mainly determined by three parameters: two principal bending constants 
and a twist constant. It is usually assumed that bending energy is isotropic. 

Recent stretching experiments [2], [3l HI |5] allow us to study mechanical 
response of a single DNA molecule. Marko and Siggia [6J reproduced the 
measured force-extension curve of DNA using the isotropic elastic rod model 
with an isotropic bending constant of about 50 nm. 

Because of DNA special structure, its bending energy is expected to be 
anisotropic. The existence of anisotropy in the bending of DNA has been 
previously reported by simulation studies as well [71 [S] . However, the exact 
values of the bending constants in the easy and hard directions (denoted 
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here by Ai and A2, respectively) are still unknown. Recently, Olson et al. 
have stated that the ratio of the hard bending constant to the easy bending 
constant is in the range of 1 to 5 [9]. 

Since the isotropic elastic rod model can explain the observed force- 
extension curve in DNA stretching experiments, one may expect that the re- 
sponse of an anisotropic DNA to the external tension is similar to an isotropic 
DNA with an effective bending constant. For a free DNA the effective bending 
constant is given by [10] 



Aeff 2 Ml A2 

We emphasize that the effective bending constant, in fact, depends on the 
external constrains applied to DNA. In case of a stretched DNA, the effective 
bending constant has been calculated by Nelson and Moroz [TT] only at the 
large force limit. In this paper, we present a perturbation theory which allows 
us to calculate the force-extension curve of an anisotropic DNA, and find the 
effective bending constant. 



2 The Model 

2.1 The Elastic Rod Model 

In the elastic rod model the DNA is represented by a continuous inextensible 
rod. The curve which passes through the rod center determines the configu- 
ration of the rod in three dimensional space. This curve is denoted by r, and 
is parameterized by the arc length parameter s (see Figure [I]). In addition, 
a local coordinate system with axes {di,d2,d^} is attached to each point of 
the rod. d^{s) is tangent to the curve r at each point 

dr* 

= (2) 

di{s) and d2{s) lie in the plane of cross section of the DNA, and are chosen 
to be in the easy and hard directions of bending, respectively. 

The orientation of the local coordinate system with respect to the lab- 
oratory coordinate system can be determined by an Euler rotation defined 

by 

A -i) = R,{-i) Ry{P) R,{a). 
a, (3, and 7 are Euler angles. The axes {di, ^2, fia} can then be related to 
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laboratory coordinate system, {x, y, z}, via equations 

di = R~^{a, (3, '-y) X , 

d2 = R-\a,P,^)y, (3) 
4 = R~^{a, (3, -f)z. 

Thus, if the Euler angles are known as a function of the arc length parameter 
s, the configuration of the rod will be uniquely determined. 
From classical mechanics we know that 

di = Qxdi i = 1,2,3. (4) 

where the dot denotes the derivative with respect to s, and Q is called the 
spatial angular velocity. The components of Q in the local coordinate system 
are denoted by Ki, K2, and uj 

= Ki di + K2 d2 + oj d^ . (5) 

These components can be expressed in terms of Euler angles and their deriva- 
tives with respect to s 



Ki = sin 7 — d sin (3 cos 7 , 

K2 = P cos 7 + a sin /3 sin 7 , (6) 
a; = 7 + d cos (3 . 

The elastic rod model introduces the elastic energy as a quadratic function 
of Q components [13] 

Eel = \kBT j ds \Ai kI + A2kI + C{uj- a;o)'] (7) 

where C is the twist constant, and ujq is the intrinsic twist of DNA. the inte- 
gral is over the entire length of the DNA. The first two terms in equation ([7]) 
correspond to the bending of DNA in the easy and hard directions, respec- 
tively. Ai and A2 are the corresponding bending constants {Ai < A2). Note 
that the bending energy is isotropic for Ai = A2. The third term indicates 
the energy needed for twisting the DNA about its central axis. 



2.2 Partition Function of a Stretched DNA 

In this section we present a standard method [HI [HI 1121 Hi] to calculate the 
statistical distribution function of the Euler angles, and to relate this distri- 
bution function to the partition function of a stretched DNA. We consider 
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here the case of pure stretching, that is, stretching with zero apphed torque. 
This situation is reahzed in many experiments [21 El IH Ej • Also we assume 
all elastic modulus sequence independent and consider them to be constant. 

Suppose that a DNA molecule is stretched by force / along z axis. Fol- 
lowing illj we neglect self-avoidance effects. Thus the DNA in our model 
behaves like a phantom chain, we also neglect the electrostatic interactions, 
which are small if the salt concentration is high enough [21 IH E] • Then, total 
energy of DNA can be written as the sum of elastic energy and the potential 
energy associated with the tensile force 



tot 



Eel - fz . 



where z is the end-to-end extension of DNA in the direction of the external 
force and is given by 

z = J -7,(13 = J COS Pds. (9) 
Using equations and 1^, one can write 

Etot = [ e{s)ds, (10) 



where e(s) is the energy per unit length of DNA and is given by 



e s 



f cos P 



'IV. 



f 

kgT- 



where / 

It is evident from equations and ffTTl) that the DNA total energy 
depends only on the Euler angles and their derivatives. This allows us to 
define a distribution function for Euler angles. For simplicity, we indicate 
the three Euler angles by the vector 9 = (a, j3, 7). In order to obtain the 
distribution function of 6, we first define the unnormalized Green function 
G(e, s I Go, 0) as follows [12j 

l-0{s)=ef 

0(9/, s| Go, 0) = / P[e]exp 
ie(o)=eo 

The path integral in ffT^ is over all paths between 0o and G/. We define 
e = T7TT , Sn = ne, and 9„ = 0(s„). Then the path integral can be written 



e s 



''ds' 



(12) 



N+l 



as 



0(s)=e 



e(o)=eo 



= lim 

e »0 

N — ►oo 
(7V+l)e=s 



Ar(e) / dQi / de 



dQ 



N 



(13) 
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where dQn = sin /5„ dan d(3n d-jn and 



Ar(e) 



^ (27re)3 J 



We call G(G, s \ Bq, 0) an unnormalized Green function since the condition 
/ G(9, s I Go, 0) dQ = 1 is not satisfied for / 7^ 0. The unnormalized Green 
function is in fact proportional to the distribution function of at point s 
for 9(0) = Go. 

The above Green function satisfies a Schrodinger-like equation [12] 



d_ 

ds 



+ H 



G(e, sIGo, 0) = 5(s) 5(0 - 60) 



(14) 



where the Hamiltonian H is given by 

•^1 



H 



+ 



/2 
•^2 



2Ai 2A 



J2 

^+iujQ J^- f cos (5 , 



(15) 



with 



J2 
J3 



—I 



—% 



cos 7 d . d ^ d ^ 

^— — + sm 7 — + cot A/ cos 7 — 
sm p oa op 07 

sin 7 d d d ~ 

^-r — + COS 7 — - cot /:/ sm 7 — 
sm/5 oa op (77 J 



(16) 



^7 



Jl, J2, and J3 are analogous to the angular momentum components of a 
quantum mechanical top with respect to a coordinate system attached to it. 
These angular momentum components satisfy the commutation relation jl5j 



[</j, Jj 



i ^ijk Jk • 



(17) 



Note that the term i uq J3 makes the Hamiltonian non- Hermit ian. In fact 
the Hamiltonian commutes with the time reversal operator and belongs to a 
class of Hamiltonians which are called pseudo-Hermitian [16j. 

The operators Ji and J2 can also be written in terms of ladder operators 

J± 

Ji = liJ+ + J-) 



J2 



2i 



:j+ - j- 
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Substituting Ji and J2 in equation ffT^ and using commutation relation ffTTj) . 
we obtain 

^ = 21"^' + ^2^ ~ 21^"^^' + Ij^-^' + + ^^0 J3 - / cos/? , (19) 
here A is the harmonic average of Ai and A2 

and 

^=4^- (21) 

A is a dimensionless parameter characterizing the anisotropy and varies be- 
tween zero and one. 

We denote the distribution function of O at the point s by \I/(B, s). From 
the definition of Green function it is obvious that "^{Q, s) can be related to 
\E'(0, 0) via equation 

^(6, s) = J G(e, sIGo, 0) ^(00, 0) dQo • (22) 

Notice that since the Green function is not normalized, \E'(B,s) is not nor- 
malized either, so we refer to it as the unnormalized distribution function. 
Considering (1^ . \I'(6,s) also satisfies equation ffT^ 

d 

HmiQ.s) = -—miQ.s) s>0. (23) 

OS 

Therefore, we can find \E'(G,s) by solving the above Schrodinger-like equa- 
tion. 

We now use Dirac notation to present our results in a more familiar form. 
Replacing \E'(6, s) with (6 |\E'(s)) we can rewrite equation (l23ll as 

^l*(^)) = -^l^(^))- (24) 

Using equations ( fT2l) . (fT3ll . and ( |22i) . the partition function of a stretched 
DNA can be written as 

Z = j {Q\^{L))dQ, (25) 

where L is the total length of DNA. Hence, in order to find the partition 
function one must solve the Schrodinger-like differential equation (l23ll and 
integrate the solution over all G values. 
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To solve equation ([221) , rewrite the Hamiltonian in the form 



where 



and 



H = Ho + XV, (26) 



^° = 21*^' + ^2^ " Ja^-^' + '^oJ3- /cos/? , (27) 



V^=^(4 + ^')- (28) 
Furthermore, we decompose Hq to its real and imaginary parts 

Ho = H^ + iH^, (29) 

where 
and 

Hl = uj^h. (31) 

ifg^ is the Hamiltonian of a quantum top. It commutes with both J3 and J^, 
where is the third component of the angular momentum operator in the 
laboratory coordinate system [T3]. Since J3 and also commute with each 
other, one can find the simultaneous eigenvectors of these three operators. 
We denote these simultaneous eigenvectors by fc, m) where k and m are 
integer numbers referring to the eigenvalues of J3 and J^, respectively. The 
quantum number n distinguishes between the eigenvectors with identical k 
and m numbers: 

H^\n,k,m)=S^^^^^\n,k,m), (32) 
J3 |n, k,m) = k \n, k, m) , (33) 
Jz \n, k,m) = m \n, k, m) . (34) 

From equations (1271) . ( l32l) . and (l33l) . it can further be seen that the eigen- 
vectors of are also eigenvectors of Hq: 

Ho |n, k, m) = \n, k, m) , (35) 

where 

^n,fc,m = ^n,A:,m + ^ • (36) 
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Since Hf^ is Hermitian, its eigenvectors form a complete orthogonal basis 
jl7j . We now expand in terms of |n, k, m) eigenvalues 

1^(5)) = e-^-^-^|n,A;,m) (37) 

and substitute |^E'(s)) into equation ([21]) • Taking the orthogonality of the 
eigenvectors into account we obtain 

g-^Cn,k,^ = -X (n,A;,m|\/|n',A;',m')e-(^"'.^'.™'-^".'=-)^C„,,fc,,™,. (38) 

n' ,k',m' 

The ladder operators in V imply that [15] 

(n, k, m\V\n', k' , m!) = {n, k, m.\V\n' , k + 2,m) Sm',m ^k',k+2 + 

(n, k, m\V\n', k-2,m) 6rn',m h',k-2 ■ (39) 

so we have 

g^Cn,k,m = -A ^(n,A;,m|\/|n',A: + 2,m)e"^^"'>'=+2.™"^°'=--^' Cn',k+2,m 

n' 

-A k, m\V\n', k-2, m)e~''^°'-^-^'^~^^^>=^"^^' Cn',k-2,m -(40) 

n' 

Substituting li'iL)) from equation ( l37|l into equation ( l25l) we can derive 
an expression for the partition function. Since J (0|n, k, m) d<d is non-zero 
only for A; = m = 0, we obtain [15j 

Z = ^4C„,o,o(i^)e-^"%'«^ (41) 

n 

where 

4 = J {e\n, 0, 0) dQ . (42) 

Thus, to determine the partition function of a stretched DNA, one needs 
to find the coefficients Cn,ofl{L) by solving the differential equation ( HOl) . 

2.3 Perturbation Theory 

In this section, we use perturbation theory to find the expansion coefficients 
and the partition function in powers of A. Let's expand Cn,k,m{s) in terms of 
A: 

00 

Cn,kAs)=Y^'^Sj'y (43) 
p=0 
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As a result, the partition function can be written as 

oo 

Z = ^APZ(f), (44) 

p=0 

where 

Z^'^ = Y.^r.C^to{L)e-'l^^^^. (45) 

By inserting Cn,k,m{.s) from equation fHHj) into equation fHUj) . one can see that 
C„^^''o(s) satisfies the following differential equations 



^^i°U = (46) 
for p = 0, and 

n' 

- ^(n, fc, m|y|n', k-2, m)e"^^"'.*-2-"^° C„^f (47) 

n' 

for p > 0. 

The value of |^E'(0)) is determined by anchoring the DNA hence independent 
of A. Thus the corresponding initial conditions are 

^iL(O) = Cn,k,m{^) for p = 

CiljO) = forp>0. (48) 
It can be seen from equation (H6l) that is constant 



n, k, m 

C^lm = Cn,k,M ■ (49) 

Therefore, the partition function to the zeroth order of A is given by 

^^°^ = EC^"'°-'' (50) 

n,k 

where 

CI = 4 Ci5,0 4,0 . (51) 

Z*^") is the partition function of an isotropic DNA with bending constant A. 
The differential equation ( H71) can be solved by iteration, and the corrections 
to Z^'^^ can be found in powers of A. The first order correction is given by 

Z^' = j:bS.e-<-^'^ (52) 

n,k 

9 



and the second order correction is given by 



Z'^' = Y.^h^:i - Ci%5,,L) (53) 



n.k 



The coefficients 6^^] and h^]. in equations ( l52l) and ( l53l) are given in appendix 
1X1 They depend on the initial conditions but do not depend on the length 
of DNA. The coefficient Url is given by 



/ J -'"1 



(m, 0|y|n2, 2) {n2, 2\V\n,0) ^ 



(CO _fO \ 

(ni, 0|\/|n2, -2) (n2, -2|y|n, 0) 



(54) 



where for simplicity, we omit the quantum number m keeping in mind that 
m = 0. Gn in equation (15^ refers to all eigenvectors with eigenvalues equal 

Clearly, if the eigenvector \n, 0, 0) is not degenerate, we have Gn = {n}. 

The coefficients bj^l are zero except for k = 0, ±2 (see appendix 1X1). 
Since the imaginary part q is kujQ, an oscillatory term with frequency 
2uq appears in Z^^\ In fact, from equation fHT]) we expect that oscillatory 
terms with frequencies {2uJo, 4a;o, ■■■ ,2puo} appear in the expression of 
Z^P\ The appearance of oscillatory terms is, in fact, an artifact of coupling 
between bending and twisting in an anisotropic DNA [18]. This is the main 
difference between the partition functions of an isotropic and an anisotropic 
DNA. Although, as we will show in the next section, this difference is not 
detectable in experiments, at least if the DNA is long enough. 



2.4 The Average End to End Extension 

Using equations (El), ( 1121) . (|22l) . and (!25l) the average end-to-end extension of 
the DNA can be calculated as [HI [H] 



<£>=i^. (56) 
L L df 

Following Marko and Siggia [B], we limit our study to the long DNA. In 
this case, because of the presence of exp(— q L) factor, the term which 
corresponds to the ground state eigenvalue of is much greater than other 
terms in the expansion of the partition function. Therefore, the partition 
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function can be approximated only by the ground state term where all other 
terms can be neglected. If we denote the difference between the ground 
state and the first excited state eigenvalues by /S.£^ then the long DNA limit 
corresponds to the condition AE^L ^ 1. We will discuss in the next section 
that this condition is indeed satisfied in the stretching experiments. 

The operator Hq is the Hamiltonian of a top in a uniform external field, 
and its ground state is unique. Thus the ground state of must be a 
simultaneous eigenvector of J3 and J2, with eigenvalues m = k = Q [6]. We 
denote the ground state and its eigenvalue by |0, 0, 0) and £^q q respectively. 
Therefore, at long DNA limit we obtain 



,(0) 



h^"' 0,0' 



and 



'^O, 0, ^ 



Since the ground state is not degenerate. Go = {0} and one can write 



(2) 



/'2 r 

^0,0,0 -'0 ; 



where 



'-'0,0,0 



(0, 0, 0|y|n2, 2, 0)1 



c-O \ 



+ 



1(0, 0, 0|\/|n2, -2, 0)1 
(^0° 0, 1 



(^o'!o,o 

Therefore, Z^^^ can be written as 

V (2) _ ^ !^ (2) _ (0) p 2 r N £ ^ 
^ — 1^^0,0 f^O.O ^0,0,0 ^ 

We expand (z) in powers of A, 



c-O N 

^n2,-2,0J 



(57) 
(58) 

(59) 
(60) 

(61) 



L L 
From equation (l56|l . we have 

L 



L 



L 



(62) 



(63) 



1 <91nZW 
^ 5/ 



(64) 



1 <9 ainz 

^9/ 



A=0 



1 a 

L df z(o) ' 



(65) 
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and 



-(2) 



) 



1 d d'^lnZ 
2L^f 



1 d 



A=0 



L 2 La/ 

Neglecting terms of order [HI [H], we obtain 



(66) 



^^0,0,0 



L 

(.(2)) 



5/ 



0, 



^^0,0,0 



(67) 



(6^ 



(69) 



L Of 

So far we have assumed that DNA is inextensible. To account for the 
extensibility of DNA the term must be added to , where B ksT is the 
stretch modulus of DNA and is about SOOkeTum"^ [6]. Thus one can write 



^^o!o,o f_ 
' df B- 



(70) 



{z ^ ') is the average end-to-end extension of an isotropic DNA with the bend- 
ing constant A. Marko and Siggia have also calculated (z^^^) [H]. Although 
they used a different Hamiltonian, i.e. 



Hi, 



2A 



/cos/3, 



our results are identical to theirs to the zeroth order. The reason is that 
and Hiso have the same ground state eigenvalues. 



3 Results 

Numerical methods are employed (see appendix [B]) to calculate the second 
order correction to the force extension curve of an isotropic DNA, assuming 
A = 50 nm, C = 100 nm [TT], and lvq = 1.8 nm~^. The result is shown in 
Figure [2j For forces slightly greater than / ~ 10nm~^ the DNA undergoes 
an over-stretching transition [12], hence the elastic rod model is not relevant. 
We have therefore picked the force range of < / < 10nm~^ to insure 
validity. 
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It can be seen from Figure [2] that (-z^^-*) is positive. Therefore, to the sec- 
ond order of A, anisotropy increases the average extension of DNA. However, 
(2;*^^-*) is small compared to {z^'^^). For A = 50 nm, the maximum value of 

the ratio ^^j^ is in the order of 10~^ (see Figure E]). 

To be sure that this result is not limited to the special case of A = 50 nm, we 
examine four different values of A in the range 5 < A < 500 nm (see Figure 
[2]). The ratio ^^joj^ for these four values of A are plotted in Figure El It can 

be seen that j^(o|| does not exceed 10~^ for A > 5nm. 

As can be seen from Figure Ej for ^4 = 50 nm, where the theoretical curve 
is best fitted to the experimental data j,6j, one must measure at least 
with the accuracy 10-^ to detect (z^^)) . Since L ~ 10 /xm in experiments 
[2], minimum accuracy of Inm is required in measuring (z). However, the 
accuracy of the experiments is by far less than this limit [2], therefore (-z*-^-*) 
can not be detected by stretching experiments. 

We now show that (z^^^) is also small. It is obvious that when \I'(B,0) 
is independent of the Euler angle 7, the partition function is invariant under 
the transformation A ^ — A. This means that odd powers of A are not 
present in the expansion of {z), i.e., (^(^p+i)) = 0. In addition, the effect of 
the initial conditions on the force extension curve of DNA is suppressed if 
DNA is long enough. As a result, one expects -^^-^^ — - to be small even when 
^(0,0) depends on 7. In other words, odd powers of A have no significant 
contribution to the end-to-end DNA extension. Therefore, to the third order 
of A, the response of an anisotropic DNA to the external tension is close to 
an isotropic DNA with the effective bending constant 

^"'=^='{hiy^ (71) 

To justify our result, we must show that the condition AS^L ^ 1 which 
corresponds to the limit of long DNA, is satisfied in experiments as well. 
Figure H] shows AS^A as a function of fA for A = 50 nm. As can be seen, 
AS^A > 1. As a result, the condition AS^L ^ 1 is equivalent to the 
condition L ^ A, which is well known in polymer physics. Since ^4 = 50 nm 
and L ~ 10 /um, this condition is satisfied in the streching experiments. 

4 Discussion and Conclusion 

It is well known that when DNA is free (i.e., no external force applied), 
the average energy of an anisotropic DNA is equal to the average energy of 
an isotropic DNA with bending constant A [8]. Moreover, Maddocks and 
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Kehrbaum [TU] have proved that in the absence of external forces or torques 
the ground state configuration of an anisotropic DNA is similar to the ground 
state configuration of an isotropic DNA with bending constant A. However, 
a stretched DNA is not free. More importantly, to calculate the average end- 
to-end extension one must deal with the free energy instead of the average 
energy or the ground state energy. 

The partition function of a stretched DNA is generally represented as 

l>00 771 

Z= v{E)exp{-j-j;)dE, (72) 

where T>{E) dE is the number of possible configuration with an energy in 
the range of E and E + dE, and Emin is the ground state energy. For an 
stretched DNA, the ground state corresponds to the configuration in which 
the DNA is fully stretched. The equilibrium configuration of the DNA is 
the configuration that minimize the free energy, F = E — ksTlu'D^E), and 
therefore is different from the ground state configuration. Clearly, bending 
anisotropy changes "D^E) for excited configurations thus changes the free en- 
ergy and equilibrium configuration of the stretched DNA. When no external 
force is applied to the DNA, the number of configurations that have the end- 
to-end extension z is exactly equal to the number of configurations with the 
end-to-end extension ~z. Consequently we have (2;) = regardless of the 
degree of anisotropy, A. Thus in the limit of fA ^ 1, anisotropy can barely 
affect the average end-to-end extension, and one expects (2;^^'') and in fact 
all the higher-order corrections to be small, as can be seen from Figures [2] 
andO On the other hand, in the limit of fA ^ 1, the energy of the ground 
state is much lower than those of the excited states, and the excited config- 
urations have a small contribution to the partition function. Therefore, the 
effect of anisotropy will be suppressed at large forces, and (-2*^^-') and all the 
higher-order corrections vanish as fA 00. This is the reason that ^^q^'^ 
is smaller at large forces (see Figure [2]) ■ 

Nelson and Moroz [H] have applied an approximate method to obtain an 
analytical expression for (z) at the limit of large forces to the second order 
of A. They found 

Aff = A{l-2{-r) = A-X'A, 
A 

with A = \{Ai + ^42), and A = |(^2 — ^i)- This result is different from 
equation (1711) . However, we rederived their calculations and obtained the 
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same result as in equation flTTj) 



Aeff = Ail - i-f) = A 
A 

Thus we believe that they just made an error in their calculations. It can 
be shown that the result of these calculations is in fact exact (see appendix 
ICj) . Therefore, at the high force limit, an anisotropic DNA behaves like an 
isotropic DNA with the bending constant A. 
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A Expansion Coefficients for Z^^^ and Z^'^^ 



Here we present the expressions for b^l and b^l- Let's use the following 
abbreviations 

and 



Vn,k,n',k' = {n, k, 0| V\n, k' , 0) , 



'-^^n,k,n',k' — '~'n,k,0 '-'n',k',0- 



The b^l coefficients are 



.(1) 



0,n,±2 ^(0) 



'n, ±2,ni,0 



, (1) _ . r K,0 ,ni,2 ^(0) 

\0 -^n 2^ 



"1 



^71, 0, ni , 2 



A^W I K,0,ni,-2 ^(0) 
"-"^1,2,0 APO '-'ni -2,0 

^^n,0,ni,-2 



and 



&S = MO, ±2. 



The coefficients are 



(,(2) 



,0,n2,±2 , ±2,n,±4 
AfO /\cO 



'-^n,±4,0i 



(73) 
(74) 

(75) 
(76) 

(77) 
(78) 
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1,(2) 

ni,ra2 



,0,n,±2 ,±2,ni,±4 ^(0) 



'-^^n, ±2,711, ±4 ^^n,±2,n2,0 



712, 0, n, ±2 ,±2,ni,0 
^<-'n,±2,ni,0 ^*-'n,±2,n2,0 



'-^ni,0,0 



(79) 



and 



f,(2) 



n,0, n2,fc ^n,k,ni,2k 



n2 fc=±2 " »ii,2fc,n2,fc n,0,ni,2fc 



^<^ni,2fc,n2,A: ^'-'n,0, ra2,fc 



^ni,2A:,0 



E E^' 



r v; 



n,0,n2,k ^n2,k,ni,0 



„j^^2fc=±2 n, 0,712, A; '"'ni, 0, 712, fc 

712, A; ^n2,k,n,0 



a 



(0) 

ni,0,0 



+ E E 



ni,n2 k=±2 
ni^Gn 



APO A CO 

^*-'7i,0,n2,fc ^'-'71,0,711,0 



^71,0,0 



E E^. 

'11,^2 
711^G„ 



r K 



n,0,n2,fc ^^2,A:,ni,0 



ni,n2 k=±2 7ii,0,n2,fe 7ii,0,n,0 



(0) 



711,0,0' 



^71,fc 



In equation 
k = -2. 





' J2k=±2 



k^O, ±2, ±4. 
indicates that one must sum over both k 



^0) 



(81) 
2 and 



B Numerical Calculations 

To calculate the eigenvectors and eigenvalues of Hf^, we use the eigenvectors 
of the angular momentum operator as the basis of the Hilbert space. We 
denote these eigenvectors by \xj,k,m)- From quantum mechanics, one knows 
that \Xj,k,m) satisfies the following eigenvalue equations \X5l 

J^\X3,k,m) = 3 U + l)\Xj,k,n.) J G Z+ U {0}, (82) 

J3\Xj,k,m) = k\xj,k,m) \k\<j, (83) 

Jz\Xj,k,m) = m\xj,k,m) \m\<j- (84) 
The vector |n, k, 0) can be expanded in terms of |xi,A;,o) as 

oo 

\n, k, 0) = '^anj,k\Xj,k,o) ■ (85) 

j=k 
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Then the equation Hf^ \n,k,0) = q |n, /CjO) transforms to the matrix 
equation 



oo 
j'=k 

where {Xj,k,o W\Xj',k,o) is given by [151 [20] 

1 . .1.1 1 



(86) 



{Xj,k,o\HQ \Xj',k,o) 



f 



v/(2j + l)(2j' + l) V(2j + l)(2j' + l) 

Similarly, (n, A;, 0|\^|n', A;', 0) can be written as 



^7) 



(n. 



k, Q\V\n\ k', 0) = J2Y1 <j,ka'n',j',k' {Xj,k,o\V\xj',k',o) , (88) 

j=k j'=k' 



where {Xj,k,o\V\Xj',k\o) is given by 



{Xj,k,o\V\Xj',k',o) 

^k,k'+2 



4A 
1 



VU + k)ij-k+ + k-l)ij-k + 2 
VU + k') U -k'+ + k'-l) (j -k' + 2) 



(89) 



The dimension of the matrix {Xj,k,o\H^\xj',k,o) is infinite. Thus, to 
solve the eigenvalue equation fISUl) numerically, we choose a cutoff for j. We 
find that the calculated values for £q ^ q and £q q q converge very rapidly. 
A choice of jmax = 120 is sufficient to calculate {z^^^) and (^'•^^) with a 
relative accuracy of 10~® (taking into account the error due to numerical 
differentiating). 



C Average End To End Extension of the DNA 
at large force Limit 

If the external tension is adequately large, the DNA remains relatively straight. 
Thus, (is lies approximately in the z direction and di and d2 will be confined. 
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as a result, in the xy plane. In this case, the components of the spatial 
angular velocity can be written in this form 



Ki = —d^x sin (f) — dsy cos (j) , 

K2 = dsx cos (f) + dsy sin , 
d(j) 



(90) 



UJ 



ds 



Where (f>{s) is the twist angle of DNA, and dsx and d^y are the components 
of (is in the x and y directions, respectively: 

ds = d3x^ + dsyf + dszZ . 

Further, since d^x and d^y are both small, we can write 



(91) 



Defining A= ^(^1 + ^2) and A= ^(^2 - Ai), the energy of DNA can be 
written as 

E = Eo + E, + Et„i,t - fL , (92) 

where 



Eq ^ 1 

knT 2 



MdL + d^y) + fidL + dly) 



ds, 



(93) 



and 



El A 



cos(20)(4^^ - d 



3y) 



ds 



A 
+ 2 



E, 



twist 



2sm{2(f)) dsxdsy ds , 



C 

— I (u — LOoY ds . 



(94) 



(95) 



On the basis of the ergodic principle, one expects that the relation 

1 r 

u{s')ds' = {00) = ujq 

holds for large s [21] • Thus, for a long DNA we can employ the approximation 
[22] 

0(s)~u;os, (96) 
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and substitute into equation to get 
El A 







A 
+ 2 



2sin(2cuos) d^xd^y 



ds 
ds . 



(97) 



To calculate the partition function, we express the total energy in terms 
of Fourier components of d^x and d^y. The Fourier transform oi d^^ + i d^y 
is given by 



j=-oo 



(9^ 



where = Using the properties of Fourier transformation, we obtain 

m 



L 



Eo = - ksTY^ [Aq] + /) 



]=~CO 



and 



A 



(99) 



(100) 



j=-ao 



where qj^ is the closest wave number to 2ujq 

o 2 7rjo 
20.0-^. 

We denote the real and imaginary parts of aj as Rj, and Ij ,respectively. 
Then the total energy of the DNA can be written in the form 



E = - fL + E^ + E^ + E, 



twist ) 



with 



E^ L °° - ~ Al °° 



'Jo- J ' 



:ior 



(102) 



j=-oo 



j=-oo 



and 



£ji L °° - ~ AL °° 

- [Aq] + f)l] + — q, qj,_j , 



fceT 2 

J = — 00 J = — OQ 

Therefore, the partition function is given by 

Z = Zr Zi ^twist ) 



(103) 



(104) 
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where 



Zn = exp(^) y H dR, exp(--^), (105) 



expl— , 



and 



-'twist 



Vlu] exp(-l^) . (107) 

The integral in equation fll07p is taken over all possible paths of uj{s). 

The average end-to-end extension of DNA can be calculated from equa- 
tion (1^ . Since Ztm^t does not depend on /, one can write 

(^) = ^-(lnZ^ + lnZ,). (108) 
df 

It is clear from equations f ll02p and (11031) that one needs to calculate only 
Zj. Zr can be calculated simply by replacing A with —A in the expression 
obtained for Zj. 

From equation (11031) we have 

E^j jo is odd. 
= { ' (109) 

E°t«i + El jo is even. 



where 

keT 2 
and 



[{M] + /) /| + {Aql_^ + /) +2AL g, g,„_, /,/,„_,] , (110) 



-^[{A + A)ql+~f]ll. (Ill) 



A;i3T 2 

For simplicity we assume that jo is odd. It can easily be shown that the final 
result does not change when jo is even. Since the variables Ij and Ij^-j only 
appear in E^^, substitution of E^ in equation (11061) yields 



fL °° r ["^ E 



:il2) 
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The integrals in equation (11121) can be calculated using the formula 



/oo 
dxdy exp[— ax^ — by"^ + 2cxy] 
-oo 



TT 



\fah 



(113) 



Then we obtain 
In 



= In Z« = ^ In Zo - i 5^ In [l - F{q,) F{q,,.,] 



(114) 



]=-oo 



where 



Zo = exp(/L) / [ J] dl.dR 



exp( 



j=-oo 



(115) 



is the partition function of an isotropic DNA with bending constant A and 



Aq' + f 



(116) 



Using equations (11081) and (11141) ,the average end-to-end extension of DNA is 
given by 



L L 2L 



^ oo 



[117) 



j=-oo 



where {z)q is the average end-to-end extension of an isotropic DNA with the 
bending constant A P, |TT] , 



(z)^ _ld\nZo 1 



L df 



:ii8) 



and 



G{q) 



A^F{q)F{q~2uJo) 



{Aq^ + f){l-A^F{q)F{q-2uo)) 
The sum in equation (I117p can be transformed into an integral as follows 

^ = ^ - ^ ljG{q) + G{2uo - q)] dq 

_ (^)o ■ 



(119) 



2n 



— oo 
oo 



Giq) dq . 



(120) 
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Changing the integration variable q to x = . — Q, and defining xq = 2. — ujq 

V / V / 

and A = 4, one can write 

r( \^ p U{x) Ujx - xq) 



with 



x^ 



1 + x 

From equations ffTTSD . (11201) and flT2Tl) we obtain 



Uix) = (122) 



T^ = 1-^=(1 + ^(A)), (123) 
2 J/A 



where g{X) is given by 



- ' ^ (x2 + 1)(1 - A2 /7(x) t/(x - Xo)) ' ^ ' 



TV 



Since Xq ^ 1 in the range of experimental data, we employ Nelson and Moroz 
approximation [TT] 

U{x-xo)-l, (125) 
to calculate the integral in equation (11241) . We find 



Thus we obtain 



g{X) = ^^=-l. (126) 



^4 = 1 ^(i-(4m . (127) 



^ 2\/fA 



A 



Comparing equation fll27p with equation flllSp . one can see that the effective 
bending constant is given by 

^.„ = A(l-(|)^)=2(-l + i-)"'. (128) 

This is the same result that we have obtained in section [H] 
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List of Figures 



Parameterization of the elastic rod. The local frame {di, d2, d^} 

is attached to the rod l26 

(color online) The contribution of the second order term to the 
relative extebtion of DNA as a function of /, for C = 100 nm, 
ujq = 1.8 nm~^. The carves show the result for different values 
of A. From top to bottom A = 5nm (Blue), A = 25 nm 
(green), A = 50 nm (black), A = 100 nm (red), and A = 

500 nm (magenta), respectively. l27 

(Color online) The ratio ^^joy^ as a function of / for the curves 
shown in Figure [21 From top to bottom A = 5nm (Blue), 
A = 25 nm (green), A = 50 nm (black), A = 100 nm (red). 



and A = 500 nm (magenta), respectively. [28 

AS^A as a function of fA for A = 50 nm, C = 100 nm, and 

c<Jo = 1.8nm-i [29 
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Figure 1: Parameterization of the elastic rod. The local frame {di,d2, 
attached to the rod. 
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f/K T (nm"') 



Figure 2: (color online) The contribution of the second order term to the 
relative extebtion of DNA as a function of /, for C = 100 nm, uq = 1.8 nm~^. 
The carves show the result for different values of A. From top to bottom 
A = 5nm (Blue), A = 25 nm (green), A = 50 nm (black), A = 100 nm (red), 
and A = 500 nm (magenta), respectively. 
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Figure 3: (Color online) The ratio |f-(o)| ^ function of / for the curves 
shown in Figure [2l From top to bottom A = 5nm (Blue), A = 25 nm 
(green), A = 50 nm (black), A = 100 nm (red), and A = 500 nm (magenta), 
respectively. 
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